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Abstract 

Phase synchronization is shown to occur between opposite cells of a ring con- 
sisting of chaotic Lorenz oscillators coupled unidirectionally through driving. 
As the coupling strength is diminished, full phase synchronization cannot be 
achieved due to random generation of phase jumps. The brownian dynamics 
underlying this process is studied in terms of a stochastic diffusion model of 
a particle in a one-dimensional medium. 



*E-mail: diego@fmmeteo.usc.es; bttp:/ / fmmeteo.usc.es 



I. INTRODUCTION 



Phase synchronization phenomena in coupled chaotic systems have been extensively stud- 
ied during the last few years in the context of non-identical chaotic systems (Rosenblum et 
al., 1996; Osipov et al., 1997), ecological systems (Earn et al., 1998; Blasius et al., 1999), 
physiological systems (Schafer et al., 1998), chaotic systems forced by an external periodic 
or noisy signal (Pikovsky et al., 1997a, 1997b), an ensemble of coupled chaotic oscillators 
(Pikovsky et al., 1997c; Osipov et al., 1997), and with an electronic model of two Rossler 
oscillators (Parlitz et al., 1996). This effect owes its name from the classical definition of 
synchronization of periodic oscillators which is described in terms of locking or entrainment 
of the phases, while the amplitudes can be quite different. Hence, synchronization of chaotic 
oscillators can be defined in the most general the locking between the phases of two 

coupled systems, while the amplitudes remain chaotically varying in time (Rosenblum et al., 
1996). 

For chaotic oscillators, there is no unique definition of phase. An approach to determine 
the amplitude A and phase of a narrow-band signal s(t) is based on the analytic signal 
concept that considers an analytical signal ip(t) as a complex function of time, ip{t) = 
s(t) +is(t) = A{t)e l<t> ^ and s(t) is the Hilbert transform of s(t) (Rosenblum et al., 1996). 
However, in other cases phase and amplitude can be defined as a function of the natural 
variables of the oscillator. For example, for the Rossler attractor <fi = arctan(y / x) (Pikovsky 
et al., 1997d) or = axctaxi{i) / ^ x 2 + y 2 ) (Rosa et al., 1998), and for the Lorenz model = 
arctan (y/x 2 + y 2 — u Q )/(z — z ) (Pikovsky et al., 1997b), where u and z are constants. 

In this paper, we focus our interest in the phenomenon of phase synchronization between 
chaotic Lorenz systems coupled unidirectionally through driving in a ring geometry. It 
has been shown, that for an appropriate set of parameters, a ring of N coupled Lorenz 
systems shows a Periodic Rotating Wave (PRW) where neighboring oscillators exhibit a 
phase difference of 2n/N and the amplitude varies with time sinusoidally (Marino et al., 
1998). This system, with a different set of parameters also exhibits Chaotic Rotating Waves 
(CRW) defined as well by a phase difference between neighboring cells of 2-k/N but the 
amplitude remains chaotic (Sanchez and Matias, 1999). In this structure there exists a 
superposition of Fourier modes k = and k = 1. Here we will show a transition from a 
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PRW with a phase difference of 2tt/N to a CRW with a phase difference of 47r/N, where 
opposite cells are phase synchronized (N even). Depending on the unidirectional coupling 
strength, random/brownian 27r-phase slips develop during the mentioned transition. 

II. MODEL 

We shall consider rings of Lorenz attractors coupled in such a way that the dynamical 
behavior (Giiemez and Matias, 1995) is defined by, 

Bj = a (yj - x j) 

yj = R (0 Xj-! + {1-P) xj) - yj - Xj Zj (1) 

4/ = x j Vj ~ b z j 

with a, R and b positive parameters. Usual parameters values are a = 10, b = | and R = 28. 
In Eq. flU), P accounts for the coupling strength, j runs from 1 to N (number of cells in the 
array), and for j = 1, x = x^- 

For p = 1, it was observed (Matias et al., 1997; Marino et al., 1998) that the synchronized 
chaotic state is stable if the size of the ring is small enough N = 2, while for a certain 
critical number N c = 3 in the case of the Lorenz model, an instability associated to the 
first Fourier mode k — 1 destroys the uniform chaotic state, leading to a PRW. As the 
size of the ring is increased, new Fourier modes become unstable and for iV = 6 a second 
instability (k = 2) develops that could lead to a Chaotic Rotating Wave (CRW) where 
neighboring oscillators exhibit a phase difference of 4n/N as it is shown in Fig. |l](a), that is, 
Fourier modes k — 1 and k = 2 compete in a nonlinear way. Thus, opposite cells are phase 
synchronized while amplitudes remain chaotic and are, in general, uncorrelated. Figure |T](b) 
shows the uncorrelated values of the amplitudes of Xj + N/ 2 (t) as a function of Xj(t). Since the 
second Fourier mode plays an important role for phase synchronization, we will focus our 
study in a ring consisiting of iV = 6 Lorenz cells described by Eq. ([]]) as P is varied. This 
phase synchronization describes the onset of long-range correlations in chaotic oscillations 
(suppression of phase diffusion), and thus also corresponds to the appearance of certain 
order inside chaos that here is shown as a CRW with certain similarities to a quasiperiodic 
motion. 
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To study phase synchronization of coupled chaotic systems, we calculate the phases of the 
oscillators and then check whether the weak locking condition A0 = \n^)j — m<f)j^.w2\ < const 
is satisfied. In this paper, we restrict ourselves to the case of m = n = 1. The definition of 
the phase for a given oscillator may be a problematic task as soon as there is not a center 
of rotation. Fig. ^| shows the (x, y) projection of an oscillator phase space for two (3 values, 
(3 = 1.0 and (3 = 0.85. At (3 = 1.0 a center of rotation can be clearly distinguished at 
(x, y) = (0, 0). Then a Poincare surface of section y = x, x > allows us to define the phase 
as (Pikovsky et al., 1997b) 

= n + -t^. (2) 

where t n is the nth crossing of the surface. Note that the phase has beeen normalized by a 
factor 27r. We see that with the surface chosen in Fig. 2(a) crossings are equal to maxima of 
the variable x(t). Therefore, we can know at what times the phase is an integer just looking 
at the time evolution of the variable x, this criterion has been used before (Blasius et al., 
1999). As long as (3 decreases unproper rotations become more frequent, (see Fig. H(b)). 
However it is clear the existence of a "rotation axis", unlike the funnel Rossler attractor 
case (see e.g. Pikovsky et al., 1997b) where an independent center of rotations emerges at 
one side of the attractor. In consequence phase is increased in one unit (i.e. 2tt radians) 
when an unproper rotation occurs. As it will be later discussed the time difference of two 
consecutive maxima does not depend on the own nature of each rotation. So it seems that 
our definition provides a "good" period. The instantaneous phase cb(t) will be determined 
through linear interpolation after calculating the instants of time at which maxima appear 
in the x(t) series. 

III. RESULTS 

The main effect of varying (3 is shown in Fig. |^. As shown above, for (3 = 1 opposite 
oscillators within the ring are phase synchronized. As (3 decreases, opposite cells still remain 
phase synchronized, except for some phase jumps. These events are defined as the non- 
occurrence of a maxima at its due time in one of the x(t) signals corresponding to cells j or 
j+N/2. In other words, we will assume that a phase jump occurs if A0(t) = <j)j(t) — <f>j+N/2(t) 
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changes by approximately ±1 for two consecutive maxima of x(t). Note for example for 
(3 = 0.76 and (3 = 0.85 two jumps are encircled in Fig. [| Unproper rotations are signaled, 
it is clear that they do not produce phase slips, although both phenomena appear when the 
signal turns more chaotic. Further decreasing the coupling strength (3 finally leads to the 
formation of a PRW with a phase difference of 1/6 (mod 1) between neighboring cells. The 
transition between PRWs and phase synchronization with jumps occurs for a critical value 
of (3 C w 0.75. 

The distribution of phase jumps is shown in the sequence of figures at the rightside of 
Fig. |H where the periods of time T between consecutive maxima of x(t) are shown for two 
opposite cells within the ring. As (3 is decreased, the map Tj +N / 2 as a function of Tj shows a 
greater dispersion from the mean value (located in the center of the figures) until the critical 
value (3 C is reached. For high values of j3, a small deviation of periods around the mean value 
appears according to the way the phase synchronization has been defined (|0i — < const 
< 1/2). As the value of (3 decreases, the dispersion around the mean value increases at the 
same time that two independent accumulation regions responsible for the phase slips appear 
(see circles at the rightside of Fig. |3|). Notice that these phase slips are not related to the 
unproper rotations which are represented by maxima (minima) peaks of the temporal serie 
that do not take a positive (negative) value (see arrows in the leftside of Fig. ||). The process 
of phase jumps formation is as follows; consecutive x(t) maxima of one cell remains phase 
synchronized with the opposite cell within the ring, until a phase jump occurs spontaneously, 
which corresponds to jumps from arms numbered (2) and (3) in Fig. ^|to the encircled zones. 
That is, phase slips are characterized by the sequences 2 — > 1 — > 2 and 3 — > 4 — > 3. At the 
same time, fluctuations in A0 (i.e. no perfect synchronization between maxima) leads to 
jumps between zones (2) and (3). Besides, it must be pointed out that as (3 is decreased 
the concept of phase synchronization defined above and used here becomes less restrictive 
as the dispersion around the mean value increases (const — > 1/2). 

The number of phase slips occurring at a given interval of time decreases as (3 is increased. 
Then, when opposite cells within the ring are phase synchronized the phase difference | A0(i) | 
is, on average, constant in time. But, if phase slips occur for /3 C < (3 < 1, then one would 
expect that for M different initial conditions, the averaged square phase difference dynamics 
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will be generally diffusive, so for large t, 

(|A0(t)| 2 ) = 2Dt (3) 

where D is the diffusion constant. Figure |] shows the linear dependence found for the root 
mean square of the phase difference (|A0(t)| 2 ) 1//2 as a function of time, in a log-log plot for 
M = 14 different random initial conditions for the Lorenz cells within the ring and for four 
different values of (3. The four graphs fit to a straight line with slope S ~ 1/2 as expected 
from Eq. ([J) (see Table I for the fitted values). 

The distribution of temporal periods r between two consecutive phase jumps is shown 
in Fig. for two different values of (3. Note the occurrence of longer periods of time r for 
higher values of (3. These distributions show an exponential decay with r as a consequence 
of the intrinsic random/brownian nature of the dynamical process underlying the formation 
of phase slips. Moreover, neither the phase jumps occur simultaneously for all couples of 
opposite cells within the array, nor the phase jumps are correlated in space, which is in 
agreement with the random dynamics of phase slip formation. 



IV. DISCUSSION 

From Fig. [| a mean value of the period (r) for each value of (3 can be defined. Now, 
by using a simple model of stochastic diffusion of a particle in a one- dimensional medium 
{random discrete walk), the averaged quadratic dispersion from the phase synchronized state 
(A0 ~ 0) is given by the following equation, 

(\Ad>(t)\ 2 ) = A (4) 

V/ 

where t/ (r) is the number of phase jumps that have appeared for t (r). Consequently, 
comparing Eqs. @ and it is possible to calculate a theoretical value for the diffusion 
coefficient D t h = (2(r)) _1 . A comparison between the diffusion coefficient D exp obtained 
after fitting the log-log plots given in Fig. [| and D t h is shown in Table I. Note the good 
agreement between both coefficients for large values of (3 as expected for a typical brownian 
dynamics. It must be noted that (r) increases dramatically with f3 (see Fig. [5] and the values 
of D t h in Table I) in such a way that it is not possible to assure the existence of an upper 
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limit of (3 above it no jumps appear. For (3 — > (3 C we have found small values of (r) of 
the order of the mean period between two consecutive maxima of x(t). Thus, jumps occur 
frequently in time and a random, uncorrelated in time, sequence cannot be assured (the 
system shows a tendence to display +1, —1, +1,... slips series). Then, the obtained values 
of the diffusion D exp are smaller than those predicted D t h using (r). 

The transition between periodic rotating waves and phase synchronized chaotic rotating 
waves has been shown to occur as the coupling strength (3 is increased. For values of /3 > j3 c , 
phase slips develop randomly in time following a diffusive process given by Eq. ([|). Note that 
the dynamics of the phase defined for a single chaotic oscillator is generally diffusive as well 
(Pikovsky et al., 1997b) and in this case, D determines the phase coherence of the chaotic 
oscillations which is inversely proportional to the width of the spectral peak of the chaotic 
attractor. On the other hand, for coupled unsynchronized nonidentical chaotic oscillators 
the average phase difference grows linearly with time (Blasius et al., 1999). Nevertheless, we 
have shown a different behavior where the root mean square of the phase difference grows 
with t 1 / 2 as a consequence of phase slips random formation. 
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FIGURES 

FIG. 1. Temporal evolution of the variable x of three contiguous Lorenz cells coupled unidi- 
rectionally through driving within a ring of N = 6 oscillators and (3 = 1 (a). The figure shows a 
Chaotic Rotating Wave (CRW) with approximate phase relationship of Att/N . Variables Xj and 
x j+N/2 (j = 1) are then phase correlated although their amplitudes remain chaotic, and uncorre- 
cted (b). 

FIG. 2. Projections onto the x — y plane of the trajectory followed by a chaotic oscillator when 
considering an array of six identical cells coupled with different values of the coupling parameter: 
(a) (3 = 1.0 and (b) [3 = 0.85. In (a) the straight line indicates the Poincare surface of section 
y = x, x > 0, while in (b) it represents the "rotation axis", since now a rotation center is not 
properly defined. 

FIG. 3. Temporal evolution of the variables x\ and x±, corresponding to two opposite cells 
within a ring of N = 6 Lorenz oscillators, (leftside) and temporal period T4 between maxima of 
x±{t) as a function of T\ (rightside) for four different values of [3. Note in the leftside of the figure 
that the phase slips are encircled. For (3 = 0.90 the arrows indicate the time intervals where the 
phase is not properly defined (here the local minima are positive). Note that these positive minima 
are not related to the encircled phase slips. At the rightside, circles mark those accumulation regions 
where the periods are greater than the average value, signalling then the occurrence of a phase 
jump. See text for an explanation of the sequence of numbers: 1^4. For [3 = 0.75 the arrow 
indicates the position of a single dot corresponding to the value of T\ = T4 = 0.36 t.u. 

FIG. 4. Log-log plot of the root mean square of the phase difference (| Ac/>| 2 ) 1//2 as a function 
of time for four different values of the coupling strength (3. Lines correspond to a linear fitting. 

FIG. 5. Distribution of temporal periods r between two phase jumps for two different values of 
(3. Note the different time scales between both graphs. Fitting lines correspond to an exponential 
decay (first order in time). 
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p 


S 


D exp x 10~ 3 (t.u.- 1 ) 


D th x 10- 3 (t.u.- 1 ) 


0.85 


0.51 ±0.02 


3.53 ±2.35 


4.41 ±0.05 


0.86 


0.47 ±0.01 


2.98 ± 1.09 


2.32 ±0.03 


0.87 


0.48 ±0.02 


1.55 ± 1.02 


0.73 ±0.02 


0.88 


0.52 ±0.01 


0.050 ±0.018 


0.055 ±0.001 


0.89 


0.48 ±0.01 


0.006 ±0.002 


0.004 ±0.001 



TABLE I: Values of the slope S, and D exp obtained from the linear fitting of Fig. |] and 
Eq. @. D th = (2(r)) _1 is calculated from the mean values of r. 
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